so (5) theory of insulating vortex cores in high- Tc materials 
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We study the fermionic states of the antiferromagnetically ordered vortex cores predicted to exist 
in the superconducting phase of the newly proposed SO (5) model of strongly correlated electrons. 
Our model calculation gives a natural explanation of the recent STM measurements on BSCCO, 
which in surprising contrast to YBCO revealed completely insulating vortex cores. 

PACS numbers; 74.20.De, 74.25.Ha, 74.60.-w 



The SO (5) model is being developed as a candidate 
theory for the high-X, superconductors. In the semi- 
nal paper by Zhangd the emphasis was on the phase 
diagram and the collective (bosonic) modes of the sys- 
tem. The first experimental support or in some sense 
the motivation for the model was the explanation of the 
41 meV magnetic resonance in the s.uperconducting state 
of YBCO observed by Mook et aln Later the energetics 
of the normal-superconductivity transition was shown ti 
explain experimental data on the condensation cnergy.l 
More recently the foundation of the SO (5) model has 
been clarified by work of Rabello et a/.cl, HenleyO, and 
Eder et aln These basic considerations provided the foun- 
dation for studying the Fermi-sector of the SO (5) model, 
and confirmed the well known electronic spectrum of the 
antiferromagnetic (AF) and the rf-wave superconducting 
(dSC) phases. However, the fermionic sector has not been 
treated in detail, and the transport properties in the nor- 
mal state still remains unresolved. 

In this paper we expand on these ideas by studying the 
fermionic sector of the S0(5) model in the superconduct- 
ing phase in the presence of Abrikosov vortices induced 
by an external Magnetic field. It has been suggested 
by Arovas et aln that in the SO (5) model the vortex 
cores can become insulating, in fact antiferromagnetic, 
in stark contrast to the standard normal metal cores of 
traditional superconductors. This remarkable prediction 
can be tested experimentally by using the improved STM 
technique to directly measure the, local electronic density 
of states in Abrikosojz. vortices.B The metal cores pre- 
dicted by Caroli et ala iiave been seen both in standard 
s-wave |S«perconductorsEl (sSC) and in the high-Tc dSC 
YBCOJiy The experimental advances naturally led to in- 
tensified theoretical studies of SC vortices. Following ini- 
tial calculations on sSIl-yortex cores ,lliHl3 the focus soon 
turned to dSC cores,E5ll3 and it was concluded that they 
are indeed metallic with states very close to the Fermi en- 
ergy. Therefore, it was a surprise when Renner et al. ob- 
served that the vortex cores of the high- Tc superconduc- 
tor BSCCO|-was completely devoid of low- lying electronic 
excitations.tll In this paper we will offer an explanation 
of this puzzling experimental observation by solving for 
the fermionic sector of the SO (5) model with insulating 
AF vortex cores. 



In our model calculation we consider strongly corre- 
lated electrons hopping on a 2D square lattice with a lat- 
tice constant of unit length. The non-interacting Hamil- 
tonian. Ha, is given by an isotropic tight-binding model: 
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points to the four nearest neighbors, and where Co-(r) 
annihilates an electron with spia, a on site r. For the 
interactions the spinor formalisnu makes it particularly 
simple to construct a SO (5) invariant Hamiltonian. In 
real space it is natural to consider the spinor 



*t(r) = {cj(r),c|(r),dT(r),di(r)}. 
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The dcr (r) operators are associated with the sites on the 
opposite sublattice of the one to which r belongs: 
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where Q — (tt, tt), and where (/3(R) is given by 
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which is only non-zero on the sublattice not including 
the origin. The long range nature of <^(R) is crucial 
for the existence of strict SO (5) symmetry. $(r) trans- 
forms like a spinor under SO (5) transformations, i.e. un- 
der rotations in the ab plane generated by the opera- 
tors Lab = iEr*^Wr"''*(r), a,b = 1,2,3,4,5, where 
P'^^ = -i [P'^,P'']. The five 4 X 4 P"-matrices are given 
in terms of tensor products of the standard 2x2 Pauli 



matrices: P^ = Go, 



'y, 



P^ = I eg) CT^, Vy ^ a,® CT„, 



T^ — 1® (Jz, and P^ — (Jx® (Jy. The indices 2, 3, and 4 
are writen as x, y, and z refering to the real spacejjlirec- 
tions of the AF order parameter. It can be shownQ that 
Li5 corresponds to the charge counting operator Q, that 
Lyz, Lzx, and L^y correspond to the spin operators S^, 
Sy, and Sz, and that Lii^x,y,z) are related to the 'ni^x,y,z)- 
operators rotating between the dSC and AF sectors. As 



in Ref. ^ we now focus on the vector interaction, which 
in the real space representation takes the form 

i/int = ^nr-s){*^(r)r'^vl/(r)}{*t(s)r«vl/(s)}. (5) 



In reahty the SO (5) symmetry is broken. However, both 
the interpretation of the 41 meV excitation as a pseudo 
Goldstone mode relating to a rotation of the dSC phase 
into the AF phase, as well as the fact that the cou- 
pling streagths in the dSC and AF sectors are almost 
identical,t2l makes it plausible that the S0(5) breaking 
is weak. A natural way to break the SO (5) is simply to 
truncate the long range correlations apparent in Eqs. (H) 
and (0) for the d-operators. 
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where (pj — (—1) ^ . Maintaining only nearest neighbors 
in the d-operator sum constitutes a simple form relating 
both-to SO (5) symmetry and to Hubbard- like models for 
dSOia. Now follow three approximations. First, we use 
the truncated d-operators instead of the d-operators in 
the interaction Hamiltonian Eq. (R). Second, we assume 
a point interaction, V{r — s) — — |V5(r — s). And third, 
we utilize the standard mean-field approximation. These 
approximations result in the following SO (5) symmetry 
broken mean-field interaction Hamiltonian: 

ifi'^f = -^y[m(r).(m(r)) + 2{A(r)(At(r))+h.c.}] (7) 

r 

with the dSC and AF order parameters given by 
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We find that in HP^l of Eq. (0) the SO (5) symmetry is 
broken in such a way that a d-wave gap function results in 
the pure SC phase, -Ek = ^k + [2V^|A|(cosfca; - cosfc^)]^, 
while a full gap develops in the pure AF phase, E^ — 



+ [^Vm{l + {cosk^ 



cosky)'^]'^. 



To elucidate the role of the gap in the AF sector we 
first study the continuum limit of our model. The im- 
portant low-lying excitations in the fermionic sector are 
concentrated in the regions near the four d-wave gap 
nodes Qa = f(cos[7r(A/2-l/4)],sin[7r(A/2-l/4)]), where 
A = 1, 2, 3, 4. We get rid of the rapid variations by 
local gauge transformations in each of the four quad- 
rants A in k-space: Ccr{r) — J^x e**^^'"?/'A,cr(i")- The gauge 
transformation is then used on Hq (with A = 0) and 
-ff j™j . Upon summing over r we keep only slowly vary- 
ing terms, i.e. terms where exp[i(QA' ± Q^) • r] van- 
ish. Not surprisingly, the only surviving terms are ei- 
ther diagonal in A or have Qa' = Qa + Q = A. This 
means that ipia{r) and '(/'3cr(r) = 'i/'icr(i") form one sub- 
space, and V'2cr(r) and ip4a{r) — ?/'2cr(r) form the other. 
Thus it becomes natural to consider the spinors ^J^ (r) = 
{V'lt(i"),V'Ax('')'V'i|(r),'0Ai (!■)}■ The gauge factor e'^^^-"" 
leads to a sign change between the terms ip\a{^ + Sj) 
and ■(/'Acr(r — ^j) in Hq and A. The difference terms aris- 
ing from this becomes derivatives in the continuum limit. 
Further care is necessary regarding extra QA-dependcnt 
signs. For simplicity we assume m(r) — m{r)ez and ob- 
tain a final Hamiltonian for the ^I^i spinor (A = 1 and 
A = 3) containing both A(r) and m(r): 
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With the Ansatz ^{(r) = (ai, 02, a^, 04)6 *'''', and k± — 
kx ± ky for k = {kx^ ky) the eigenvalue problem becomes 



0, (11) 



A pure SC phase has a constant SC order parame- 
ter A, while 771 = 0, and the spectrum becomes E = 

■iiJt'^k\ + |Apfc^. The corresponding eigenstates are 



easily found. For a pure AF phase tti is a constant and 
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A = 0. The spectrum now becomes E = ±1 /t-^fcl + m^, 
with associated eigenstates. 

We now imagine the plane to be divided into two parts. 
For X < the system is in the SC phase while for x > it 
is in the AF phase. It is now simple to study the scatter- 
ing problem where a particle with energy |_E| < 777 in the 
SC sector is moving towards the barrier formed by the 
AF sector. The result is not surprising: if the particle 
starts out with a momentum near, say, Qa it is reflected 
completely by the AF sector (where it only acquires an 



exponentially damped probability) , and it ends up with a 
momentum near either Q^ or Q^. The process resembles 
Andreev reflection in the quantum number A. The con- 
clusion of this exactly solvable model is clear: low energy 
particles in the SC sector can be confined by a surround- 
ing AF sector, or conversely, the AF sector expels low 
energy particles. 

We now proceed to discuss dSC vortices, first briefly 
mentioning the case of normal cores followed by our 
S0(5) model calculation of vortices with AF cores. In his 
semiclassical analysis of the electronic densite of states 
produced by d-wave vortices, Volovik showedllj that only 
a small part of the density of states results from quasi- 
particles localized at the vortex cores, and that that 
part is a function of the vortex density. Hence, in any 
realistic calculation of quasiparticle states, the entire 
vortex lattijee must be taken into account. Wang and 
MacDonaldM made the first self-consistent, numerical 
lattice calculation of a tight-binding model for d-wave 
type-II superconductors using the Bogoliubov-de Gennes 
equations. In the following we expand their work to the 
case of d-wave superconductors with antiferromagnetic 
cores as described by the SO (5) model. We take ad- 
vantage of earlier self-consista*Ht calculations of isolated 
vortices in the S0(5) model,0Jl3 where in the symmetric 
gauge Ao(r) = A °'°^^> eg the SC and AF order param- 
eters as a function of the distance to the vortex core in 
polar coordinates are given by A(r) = /o(?')e*^*-''-' and 
m{r) = TOo(r), respectively. Here the functions foif)i 
ao{r) and mo{r) are only known numerically. Note that 
in(r) points in a constant direction, so only the size 
m{r) is given. For a lattice of non-overlapping vortices, 
i.e. vortices further than a few times the London length 
apart, the self-consistent solutions for vector potential, 
the SC and the AF order parameter are expressed in 
term of the vortex centers, Rj, as 

A(r) = /(r)e^^M ^Hfoir - R,)eE.^Arg(r-R.) ^^^^ 

j 
m(r)=^mo(r-R,) (13) 

j 
A(r)=^Ao(r-R,), (14) 
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where Arg(r — R^) is the polar angle between r and Rfc. 
In a lattice model a particularly simple way to construct 
the magnetic unit cell is the following. For each area 
penetrated by one flux quantum h/e a_pirac anti-vortex 
string carrying a flux —h/e is added.tJ The strings will 
have no physical consequences at all when placed between 
lattice sites. However, they allow for the construction 
of a vector potential periodic in the magnetic unit cell, 
through which the magnetic flux is zero. 

We now construct a square Abrikosov lattice with a 
Dirac anti-vortex added to the center of every second 
vortex. Since each vortex carry half a flux quantum, the 
smallest magnetic unit cell possible consists of two vor- 



tices. However, due to better convergence properties in 
obtaining the periodic vector potential and a periodic 
representation of the SC order parameter (especially its 
phase 0{v)), we choose to double the magnetic unit cell. 
Our unit cell contains two ordinary vortices on one diag- 
onal and two vortices penetrated by Dirac anti-vortices 
on the other. Periodic forms of A(r), 0(r), and m(r) 
are then easily found by adding up contributions from a 
large number of unit cells (typically 64) surrounding the 
one we are studying. From this we obtain a mean field 
lattice Hamiltonian H = Hq + H^l given by Eqs. dll) 
and Im), which is periodic in our unit cell. Based on the 
Bogoliubov transformation for operators within our unit 
cell 



(7") 



a\t — 



{u^ir)cUr) + av^ir)c,{r)}, (15) 



where ct = ±1 is the spin index and a — —a, the equa- 
tion of motion for the 7"-operators using the periodic 
Hamiltonian H leads to the Bogoliubov-de Gennes equa- 
tion for the eigencnergies and eigenstates of the fcrmionic 
quasiparticles: 



T + aM D 

D* -T* + aM 



E° 



(16) 



Here E"^ is the quasiparticle energy, u" and v" are vec- 
tors containing the values of M"(r) and w"(r) on each 
lattice site in our unit cell, while the block matrices T, 
£), and M are given by 



(T)rr'=-te-^^/>'^'^<5r',r+5, 



J = l 



(D)rr'=^^,[D(r')+I?(r)](5r',r+5, 
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(M)rr' = ^ ^jV3j'M(r'-5j)(5r',r+5,+5^, 

-M(r')^r',r, (19) 

with D{y) = iF(A(r)} and M(r) = ie*'^"'V'(m(r)). 

In the numerical calculation we use a N x N lattice 
with N — 44. The origin is put in the center and the four 
vortices in the center of each of the quadrants. The peri- 
odicity is ensured by having H{r + Ndi) = H{r+N62) = 
H{r). The Bogoliubov-de Gennes equation, Eq. (16) be- 
comes a 2N'^ X 27V^ eigenvalue problem yielding for a 
given value of the spin variable a the spectrum E"' and 
the Bogoliubov coefficients u" and v". To compare our 
calculations pwith the experimental STM measureipents 
on vorticeg3^t30 and with the existing calculationsOJ on 
ordinary sSC and dSC vortices we compute the tempera- 
ture dependent local density of[States-|(LDOS) according 
to the standard minimal modelllilO'til 



a. 

+ |«"(r)|2{-/'(£;"+i?)}], (20) 

where /(e) — [exp(e/fcBT) + 1]^"'^, and where we have 
neglected the dispersion in the magnetic BriUouin zone. 
The calculation yields the LDOS shown in Fig. ^ In all 
cases V = 0.8i, k-^T = O.li and ji = —0.6t, which due to 
the band structure leads to an asymmetric LDOS. 
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FIG. 1. LDOS in the vortex core (full line) and in the bulk 
SC (dashed line) for (a) BCS s-wave SC, (b) BCS d-wave SC, 
and (c) SO (5) SC with an AF vortex core. 



half filling, in contrast to the hole doped bulk material 
maintained at lower filling. Such a charging energy must 
be taken into account in a detailed calculation of the 
energy gained by forming an AF vortex core. Our cal- 
culation of the generic features in the fermionic sector 
of the SO (5) model shows that the measured LDOS can 
be explained if one simply assumes that YBCO with its 
metallic vortex cores is a pure dSC SO (5) superconduc- 
tor, while BSCCO is a dSC/AF SO (5) superconductor. 
We obtained our results by studying both the analytically 
solvable model of a perfect SC/AF interface and by exact 
numerical diagonalization of an Abrikosov lattice model. 
Clearly, further theoretical insight in the dual dSC/AF 
nature of the high- Tc compounds can be obtained from 
studies of the striped phases, whewj alternating stripes of 
SC phases and AF phases occurE2l. 
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First, to check our calculations, we change the model 
from S0(5) to ordinary sSC and dSC. The latter is pro- 
duced by setting Mr,r' = in Eq. ( p^ , and the former by 
furthermore setting (-D)r.r' — ^V{A{r)}Sr,r' in Eq- (18). 
As shown in Fig. nk and n b we confirm qualitatively the 
main conclusions of Refs. 14|-[l6|. In the bulk of the sSC 
phase a full gap is observed, while a mid-gap peak (which 
splits at T = 0) develops in the center of a sSC vortex. 
In the bulk of the dSC phase a steady rise of the LDOS is 
seen around the mid-gap position, while a mid-gap peak 
develops in the center of a dSC vortex. Our model calcu- 
lation captures mainly generic features and can therefore 
not be used in the ongoing debate Qt|t|ie detailed form 
of the LDOS in the dSC vortex coreJialla However, this 
issue is not important for our main observation in the 
SO (5) case: instead of a mid-gap peak the LDOS is dra- 
matically suppressed in the AF vortex core resembling 
bulk behavior as shown in Fig. ^c. This confirms the 
conclusion of the dSC/AF interface in the SO (5) model 
studied in the first part of this paper. The AF phase 
effectively suppresses any fermionic low energy states. 

We thus reach our main conclusion. The experimen- 
tally observed lack of electronic quasiparticle states in 
the center of Abrikosov vortices in BSCC0ll3 as opposed 
to the measurements of a normal metallic core of vortices 
in YBC0l2I finds a natural explanation in the framework 
of the SO (5) model. As already pointed out by Arovas et 
al.El the nature of the SO (5) vortex cores are governed by 
the parameters {e.g. doping level and coupling strengths) 
of the given high- Tc material. The cores can either be- 
come metallic, i.e. a pure dSC behavior, or insulating, 
i.e. a mixed dSC/AF behavior. At the present stage of 
the S0(5) theory it is difficult to predict which materi- 
als will in fact develop AF vortex-jCores. For example, 
as is studied in the striped phaseE3, the insulating vor- 
tex cores are negatively charged, since they must be at 
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